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We study the temperature dependence of optical properties of dilute gases of alkali-metal atoms 
in the state with Bose-Einstein condensates. The description is constructed in the framework of 
the microscopic approach that is based on the Green-functions formalism. We find the expressions 
for the scalar Green functions describing a linear response of a condensed gas to a weak external 
electromagnetic field (laser). It is shown that these functions depend on the temperature, other 
physical properties of a system, and on the frequency detuning of a laser. We compare the relative 
contributions of the condensate and non-condensate particles in the system response. The influence 
of the temperature effects is studied by the example of two- and three-level systems. We show that 
in these cases, which are most commonly realized in the present experiments, the group velocity 
and the absorption rate of pulses practically do not depend on the gas temperature in the region 
from the absolute zero to the critical temperature. We discuss also the cases when the temperature 
effects can play a significant role in the phenomenon of slowing of electromagnetic pulses in a gas 
of alkali-metal atoms with Bose-Einstein condensates. 
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I. INTRODUCTION 

Bosc-Einstcin condensate (BEC) is one of the most 
impressive examples when the matter demonstrates its 
quantum nature on the macroscopic level. Now this sys- 
tem is interesting also due the possibility of observing the 
electromagnetic pulses propagating with extremely slow 
group velocities in it [l[. 

Up to now, in the theoretical investigations pretending 
to describe the mentioned phenomenon in a BEC (see 
Refs. @, [|[) the authors assumed that the temperature 
of a gas is small in comparison to the critical tempera- 
ture. In other words, this assumption corresponds to the 
consideration of the zero-temperature limit. But in real 
systems, the temperatures can be of the same order as 
the critical temperature. Therefore, one needs to study 
the account of temperature effects in the ultraslow light 
phenomenon. Naturally, it is also important to compare 
the theoretical results to the experimental data [l[ de- 
scribing the dependence of the group velocity of a signal 
on the temperature of a system. 

In the present paper, we generalize the approach devel- 
oped earlier in Ref. 0] for the uniform (nontrap) systems 
on the case of finite (nonzero) temperatures. This ap- 
proach is based on the Green-functions formalism Q and 
an approximate formulation of the second-quantization 
method [5]. An object of the mentioned generalization 
is a study of the influence of temperature effects on the 
dispersion characteristics of the system and, as a result, 
on the propagation properties of a signal in it. 



II. LINEAR RESPONSE OF A GAS IN A BEC 
STATE AT FINITE TEMPERATURES: GREEN 
FUNCTIONS 

To describe the optical properties of gases consisting 
of alkali-metal atoms that are used in the BEC-related 
experiments, it is most convenient to use the model of an 
ideal gas of hydrogenlike atoms in the stationary state 
(the limits of this approach are discussed in Ref. (H). 
By the term "stationary state" we mean that the atoms 
can be found only in the states whose lifetimes are much 
greater than the relaxation time of the system (e.g., hy- 
pcrfinc levels of the ground state) and in the states whose 
occupation are stimulated by an external electromagnetic 
field (e.g., a laser radiation). In the case of a BEC pres- 
ence in this gas, the density distribution of atoms in the 
quantum state a by the momentum p at nonzero tem- 
peratures (0 < T < T c , T is the temperature in energy 
units) can be set equilibrium, therefore, it can be written 
as follows (see also Ref. [4]): 



fa(p) [l - (T/T CQ ) 3/2 j 5(p) 

+ 5a (2^)- 3 {exp [(e a (p) - M a)/T] - 
where v a is the total density of atoms in the a state 

dpv a (p), 



(1) 
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where T ca is the temperature of the transition of a gas of 
atoms in the a state to the BEC phase, 5(p) is the Dirac 
delta function, and g a — (2F a + 1) is the degeneracy 
order of the state by the total momentum F a of an atom 
in this state. In the next calculations, we assume that an 



2 



external field is present in the system, i.e., we set g a — 1. 
Here, also h is the Planck constant, e a (p) = e a +p 2 /2m, 
where e a is the energy of an atom in the a state, m is 
the atomic mass (m = m p + m c \ m p and m c are masses 
of the atomic core and electron, respectively), and fi a is 
the chemical potential an atom in this state. 

Now we must note the following. It is known that 
a Bose-Einstein condensation is a collective effect. The 
atoms occupying the same quantum state a are identical. 
Because of this fact, one cannot definitely state what par- 
ticles participate in the condensation process. Therefore, 
the terms "condensate" and "noncondensate" particles 
do not correspond directly to the selected atoms in the 
system. Strictly speaking, the term condensate particles 
corresponds only to the fraction of all atoms, which par- 
ticipate in the formation of the Bose-condensed compo- 
nent. Hence, the number of noncondensate particles can 
be found from the difference between the total number of 
atoms and number of condensate particles in the system. 
By the use of these terms, now we can say that the first 
summand in Eq. ([T]) corresponds to the contribution of 
condensate particles in the distribution function and the 
second summand corresponds to the normal component. 
Thus in the next calculations, we differentiate contribu- 
tions from the different types of particles in the response 
of the system to the external perturbation by an electro- 
magnetic field. 

It is shown in Ref. [|| that the linear response of an 
ideal gas of hydrogcnlike atoms with BECs to a weak 
electromagnetic field can be studied from the first prin- 
ciples in the framework of the Green-function formalism. 
There it was shown that the Fourier transform for the 
scalar Green function can be defined by the relation (see 
also Ref. @ for details) 



G(k,c) = i]T5> Q/5 (k)| 2 

P a,P 
/q( P -k)-/g(p) 

£ Q (p) - Sp(p - k) + U + ij al 3 ' 



(2) 



Here, by k and u> we denote the wave vector and the 
frequency of the external perturbing field, respectively, 
and -f a p denotes the linewidth related to the probability 
of a spontaneous transition between the (3 and a states. 

The quantity <j a p (k) defines the matrix elements of the 
charge density of hydrogenlike atoms. This quantity can 
be expressed in terms of the wave functions ip a and jpp 
of the atoms in a and (3 states (see also Refs. [3, [5|), 
respectively, 



CTa/3(k) 



' / rfy^(y)^(y) 

(^ky 



exp 



)- eXP ("V ky ) 



where e is the elementary charge. In particular, in the 
case when the dipole transition is allowed between the a 
and (3 states, in the linear order over small term ky <C 1, 



one gets 

CT Q/3 (k) ikd Q/3 , d a0 = e J dyy(p^(y)(pp(y), 

where d Q ^ is the matrix element of the atomic dipole 
moment. It should be mentioned that in next calcula- 
tions, we use only the quantities that are proportional to 
|cr Q ^(k)| 2 . Therefore, for the allowed dipole transitions, 
these quantities can be expressed in terms of the average 
dipole moment d a /3, 



\a af3 (k)\ 2 a k 2 d 2 af3 /3. 



(3) 



The quantity / Q (p) in Eq. corresponds to the Bosc 
distribution function of atoms by the momentum p, 



f a (p) = {exp[(e Q (p) - fi a )/T] - 1}- 



(4) 



Note that one can use the Bose distribution in the 
form Q only in the case of the thermal equilibrium in 
the system. Thus, in the next calculations, we consider 
only the states whose lifetimes are much greater than the 
system relaxation time and the states whose existence in 
the system is stimulated by an external field. Hence, 
for a studied gas (a condensed phase can be formed by 
atoms with the energy e Q ), one gets the condition for the 
chemical potential fi a (see Ref. |7| for details) 

fJ> a (T < T c ) = e a . 

At temperatures < T < T c in Eq. ©, one can sub- 
stitute the distribution function ([4]) by the density dis- 
tribution |T|). To this end, one needs also use the rule 



Next let us note that the relation ((2|) is linear over the dis- 
tribution function / Q (p). Therefore, in accordance with 
Eq. {1} , it can be divided into two summands that define 
the contributions of the condensate and noncondensate 
particles in the system response 



G(k, u) = G (c) (k, u) + G (n) (k, w), 



(5) 



where the summands G^ c ^(k, uj) and GW(k,w) are de- 
fined by 



GW(k, W ) = 2|^(k)| 



a,/3 



2 ( 1 ^ck ) 



G (n >(k, w ) = ^r^i^k)! 2 

a,f3 

oo 1 

2itp 2 dp f dy 



(6) 



exp (e p /T) - I J Suj a fj + pky/m + ij a/ 3 ' 

-1 

Here, t a — T/T ca denotes the relative temperature of 
a gas (in this paper, we consider the case t a < 1), 
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8oj, 



a/3 



uj + Ae a/ 3 is the frequency detuning taken in the 

„2 , 



energy units, Ae Q( g — e a — sp, e p — p /2m denotes the 
kinetic energy of an atom, and the integration variable y 
denotes the cosine of the polar angle 9, y = cos8. Here 
and below, we neglect of the recoil energy e r = h 2 k 2 /2m 
that is small enough to do a significant contribution into 
the effect (e r <C 7q /3; strictly speaking, it can be ac- 
counted by redefining the quantity 8ui a p, i.e., by shifting 
the resonant frequency by e t ). For the definiteness, in 
the next description, the summand G^ c ^(k, u>) in Eq. ([5]) 
is called as the condensate Green function and the sum- 
mand G< n )(k, uj) is called as the noncondensate Green 
function. 



A. Condensate Green function at finite 
temperatures 

Firstly, we consider the scalar Green function corre- 
sponding to the contribution of condensate particles. For 
a convenience, let us study in detail the real and imagi- 
nary parts of it. It is shown below that in the region of 
transparency the real part makes a main contribution to 
the refractive index of a gas and imaginary part makes a 
main contribution to the absorption rate of light pulses. 
In accordance with Eq. ((6]), the real and imaginary parts 
of this function can be written as 



ReG< c > = £)Ks(k)| s 
ImG< c > =-£> a/3 (k)| : 



((W Q/3 ) 2 + (7a/3) 2 ' 

(<5w Q/3 ) 2 + (7c/3) 2 ' 



(7) 



To simplify the description, in the next calculations, we 
use a model of a two- level system (a = 1, j3 = 2). Below, 
by the index "1" we denote the set of quantum numbers 
corresponding to the ground state and by the index "2" 
we denote the set of quantum numbers corresponding to 
the excited state. It is also assumed that the occupation 
of the exited state is stimulated by a low-intensity laser 
pulse. Due to a low intensity of the pumping field, we 
can consider that the density of the atoms in the excited 
state is small in comparison to the density of atoms in 
the ground state (see Ref. @ for details) 

V2 < v\ = v. 

According to this relation in Eq. ([7]) we neglect of 
the summands that are proportional to the density v^,. 
Therefore, one can get 



ReG<°> = |a 12 (k)| 2 Kl - t 3 ' 2 )^ 1 F' {A), 
ImG^ = |a 12 (k)| 2 z/(l - i 3 / 2 ) 7 - 1 F"(A), 



(8) 



where we introduce the quantity 7 = 712 and the di- 
mensionless functions F'(x) = x/(x 2 + 1) and F"(x) = 
— l/(x 2 + 1). The dependencies of these functions on the 
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FIG. 1: (Color online) Dependencies of the functions F' and 
F" corresponding to the response of the condensate particles 
on the frequency detuning A. 



relative frequency detuning A (A = Slo/j) are shown in 

Fig.m 

Let us note that the functions F' and F" do not de- 
pend on the temperature. Hence, it is easy to see from 
Eq. ijHJl that the real and imaginary parts decrease with 
temperature by a power law as t 3 / 2 . Note also that the 
condensate Green function is proportional to the density 
of condensate particles in a gas, = 2/(1 — t 3 / 2 ), and 
it equals to zero at the critical temperature, t = 1. 



B. Noncondensate Green function at finite 
temperatures 



Now we study the scalar Green function corre- 
sponding to the non-condensate particles (see Eq. ©). 
Taking the real and imaginary parts, for the two-level 
system (setting the indexes a — 1 and j3 = 2, see above), 
we come to the relations 



ReG^ = |a 12 (k)p 
ImG (n) = |a 12 (k)| 2 



, m 2 T 
Aw 2 h 4 k 
m 2 T 

4ir 2 h 4 k 



/'(A,t), 
/"(A,t), 



(9) 



where we introduce the dimensionless functions F(A,t) 
and F'(A,t) that depend on the frequency of an exter- 
nal field and on the system parameters (including the 
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temperature of a gas) 

oo 

l'(A,t) = f x ' h ' 



- 1 



In 



1 + [A + x{nVi)}' 



I"(A,t) = 2 



xdx 



-A 



1 + [A-x(nyft)} 2 
arctan A — x(nVi) 



(10) 



arctan 



A + :r 



(kV*)] } 



Here, analogously to the formulas for the condensate par- 
ticles (see Eq. (|SJ)), A = 5ui/j is the relative frequency 
detuning, t = T/T c is the relative temperature, 



T c = 2 7 r[C(3/2)]- 2 / 3 ftV/ 3 m- 1 



(11) 



is the critical temperature of the transition to the BEC 
phase, and C( x ) is t ne Riemann zeta function. In 
Eq. (|10| . we also introduce the dimensionless parame- 
ter K, 



hk l2T c hk ,_ i ., 

k = — \/ — hv 1 ' 3 

7m 



(12) 



that depends on the physical properties of a gas and the 
wave number k of the external electromagnetic field. It 
is shown below that the value of this parameter defines 
the characteristics of the normal component response to 
an external perturbation. 

Let us estimate the value of the parameter k for the 
resonant radiation that corresponds to the sodium D2 
line. The laser pulses tuned to the components of this 
line were used in Ref. [l[. Taking m — 3.82 x 10~ 23 g, 
7 = 8.1 x lCT 21 erg, k = 1.07 x 10 5 cm" 1 , and T c = 
435 nK, we get k ~ 0.025. Using this parameter value it is 
easy to get the dependencies for the integrals I' (A, t) and 
I"(A,t) that are shown in Fig. O There, one can study 
also a behavior of the integral I' with a decrease of the 
temperature to the values 0.5T C and 0.1T C , respectively 
(note that the integral I" has an analogous temperature 
dependence) . 

From Fig. [21 one can see that the dependencies of the 
functions I' and I" are similar to the dependencies of the 
condensate functions F' and F" (cf. Fig. [TJ). The rea- 
sons for this similarity become evident below, but now it 
is easy to see that the values of the functions V and I" 
strongly depend on the gas temperature. In particular, 
from Eq. © and Fig. [21 we can conclude that the con- 
tribution of the non-condensate Green functions in the 
total Green function ([5]) increases with temperature. Let 
us study this effect in detail in the next subsection. 
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FIG. 2: (Color online) Dependencies of the functions I'(A,t) 
and I" (A, i) on the frequency detuning at T = T c (bold lines). 
By thin lines, the dependencies of the function I' at temper- 
atures T/T c = 0.5 and 0.1 are shown (k = 0.025). 



total response of a gas in the case k«1 can be simplified 
at arbitrary values of the detuning A. Really, due to the 
strongly decreasing function exp(— x 2 ) in the integrand, 
the main contribution to the integral give small values of 
the integration variable x, x < 1. This fact allows us to 
expand the integrands into series over (xk^/I) <C 1. As a 
result, using Eqs. ([5]) and (|10p , accurate within quadratic 
summands, we get 



-F' 2 (A) + F"(A) 
3 



J'(A,t)«W*F'(A)j|-^ 
+(kV*) 

I"(A,t)^A K VtF"(A)y^- i 
According to the relations 



(13) 



F ,2 (A) - V' 2 (A) 



x dx 
e x2 - 1 



x dx 



1 4 



-C(3/2), 



x dx 3-^7? 



C(5/2), 



C. Dependence of the scalar Green function on the 
temperature 

Let us note that the functions I' and I" characterizing 
the contribution of the non-condensate particles to the 



the Green functions (0 with account of Eqs. [) — (p" 
can be written as 



ReG^ w |ai 2 | 2 ^ 3 / 2 7- 1 F 1, (A)[l + K 2 t5'(A)], 
ImG (n) w | ( ti 2 | 2 ^ 3/2 7- 1 F"(A)[1 + K 2 ^"(A)], 



(14) 
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FIG. 3: (Color online) Relative dependencies of the conden- 
sate and noncondensate Green functions on the temperature. 
All functions are normalized to a value of the condensate 
Green function at t = and corresponding detuning A. 
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FIG. 4: (Color online) Relative dependencies of the total 
Green functions on the temperature. All functions are nor- 
malized to a value of the condensate Green function at t = 
and corresponding detuning A. 



where we introduced the functions 



S'(A) 



, 05/2) 
J C(3/2) 



_ 3 05/2) 
S ( } ~ 2 C (3/2) 



F' 2 (A) + -F"(A) 



F' 2 (A)- -F" 2 (A) 
3 



Note that Eq. (fl4|) explains also the similarity in the 
behavior of the functions F' , F" and I" at (n\ft) < 1 
that is mentioned above (see Figs. Q] and [2]). 

It is easy to see that the "total" Green function in the 
region from the absolute zero to T c weakly depends on 
the temperature. This dependence appears only in the 
quadratic terms over {n^/t) <C 1. Really, by the use of 
Eqs. ©, ©, and HID, one can get 



ReGn\<T 12 \ 2 vj- l F'{A)[l + K 2 tS'{A)}, 
ImG« la^^F'^A^l + K 2 tS"{A)}. 



(15) 



But in the case (nVi) > 1, the expansions (fl3|) become 
incorrect and analytical formulas for the functions /' and 
I" cannot be found. Therefore, in this case, we use the 
numerical calculations at the definite values of the de- 
tuning A and parameter k (see Figs. [3J and [4}. 

In Fig. [3J one can see the dependencies of the conden- 
sate and noncondensate Green functions on the tempera- 
ture. The real parts of these functions are normalized to 
a value of the real part of the condensate Green function 
at t = and corresponding detuning A. Analogously, 
the imaginary parts of these functions are normalized to 
a value of the imaginary part of the condensate Green 
function at t = and corresponding detuning A. In 
other words, in Fig. [3J the following dependencies are 



shown: 

Re G (c ' n) (t, A) / Re G (c) (t = 0, A), 
Im G (c ' n) (t, A) / Im G (c) (t = 0, A). 

One can see that due to this normalization, the inequality 
Im G(t, A) / Im G (c ) (0, A) > takes place. But note that 
the imaginary part of all Green functions is negative due 
to the function F"(A) (see Eqs. © and (fT4l» . 

Analyzing the dependencies shown in Fig. [3J one can 
conclude that the real and imaginary parts of the con- 
densate Green function demonstrate a similar behavior 
(drop-down bold curve). This dependence agrees with 
the obtained relations. But at the same time, the rela- 
tive values of the real and imaginary parts of the non- 
condensate Green functions strongly depend on the pa- 
rameter k and detuning A. One can conclude that at 
k < 1 , by the use of the mentioned normalization, the 
real and imaginary parts are practically coincide. In the 
opposite case, k > 1, as one can see in Fig. [3j the real 
and imaginary parts demonstrate different dependencies. 
This fact must have a strong impact on the dependence of 
the total Green function (the sum of the condensate and 
non-condensate Green functions) on the temperature. 

In Fig. HI the dependencies of the total Green function 
on the temperature are shown at different values of the 
parameters A and k. Here, it is important to pay at- 
tention to the central curve that is practically horizontal 
(the case k< 1). This curve demonstrates a weak de- 
pendence of the total Green function (including its real 
and imaginary parts) on the temperature. Note that one 
can come to the same conclusion by analyzing Eq. (|15p . 
Thus at k <C 1, the influence of the temperature effects 
in the response of a BEC to an external field is insignif- 
icant. The problem on obtaining the response of this 
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system can be solved with a good accuracy by setting 
the temperature of a gas equal to zero (exactly this ap- 
proximation was used in Refs. 0, IE H, [§] ) ■ In other cases 
(at k > 1), as one can see from Fig. @J the influence of 
the temperature effects can be rather significant. 

Evidently, the effects studied in this section in some 
way (significantly or not) can influence on the dispersion 
properties of a condensed gas. Let us recall that the res- 
onant peculiarities of the dispersion characteristics play 
an essential role in the ultraslow-light phenomenon in a 
BEC of alkali-metal atoms. 



III. DISPERSION CHARACTERISTICS OF A 
GAS IN A BEC STATE AT FINITE 
TEMPERATURES 
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It is known that in the framework of the linear ap- 
proach, the permittivity e(k, ui) of a gas can be expressed 
in terms of the scalar Green function (see in this case 
Refs. HH) as 

e- 1 (k, W ) = I + ^G(k,c). 

Therefore, the real and imaginary parts of the permittiv- 
ity e(k,w), 

e(k,w)=e / (k,w)+*e // (k,aO 
can be written as 

1 + 4vrfc- 2 ReG 



6 ^' ^ (1 + 47rfc- 2 Re G) 2 + (47rfc- 2 Im G) 2 ' 
„ -47rfc~ 2 ImG 

6 ( ,UJ ' = (1 + Airk- 2 Re G) 2 + (4irk~ 2 Im G) 2 ' 



(16) 



At the same time, the refractive index and damping fac- 
tor can be written in terms of the real and imaginary 
parts of the permittivity, 



l'(k,w) = -^\/Ve' 2 +^ 2 + e', 



(17) 



Hence, basing on the derived equations for the scalar 
Green functions, one can study the propagation prop- 
erties of weak electromagnetic pulses through a BEC at 
finite (nonzero) temperatures. As the derived analytical 
expressions for the refractive index and damping factor 
in the general case may have rather lengthy form (see 
Eq. ((S])), let us consider some particular cases. 



A. Temperature effects in two-level systems 

Firstly, let us consider a case when the frequency of an 
external field is close to the energy spacing between two 



FIG. 5: (Color online) Dispersion dependencies of the two- 
level system in the BEC state. The used value for the param- 
eter a = 0.015. 



definite quantum states of atoms. Thus, using the most 
common method (see, e.g., Ref. in calculations we 

can consider only the resonant terms corresponding to 
this transition. In other words, we can consider atoms 
as a two- level system. Therefore, by the use of Eqs. ([5]), 
©, ([!]), and (jT2J) , the real and imaginary parts of the 
permittivity (see Eq. (I16p ) take the form 



e'(k,cj;i) = 
e"(k,«;t) = 



1 + aGx 



(l + aG 1 ) 2 + (oG 2 ) 2 ' 

-aG 2 
(l + aG 1 ) 2 + (aG 2 ) 2 ' 



(18) 



where 



d = (1 -t 3 / 2 )F'(A) + btI'(A,t;b), 
G 2 = (l-t 3 ^ 2 )F"(A) + btl"(A,t;b), 



47T 



ki 2 (k)| 2 ^ 
fc 2 7 



(19) 



(20) 



In particular, for the condensed sodium vapor with the 
density v — 1.44 x 10 12 cm~ 3 that interacts with the 
resonant radiation corresponding to the Z) 2 line, in ac- 
cordance with Eqs. Q and (|20)) . we get a w 4.19d 2 z//7. 
Next, considering d 2 w Sff> (3.52er ) 2 (r is the Bohr ra- 
dius, e is the elementary charge, and Sff 1 is the relative 
strength of the dipole-allowed transition F — > F' [H|), 
F = 2, F 1 = 2, and S" 22 = 1/4, we come to a w 0.015 and 
b w 20.38. For the two- level system with these parame- 
ters, the characteristic dependencies are shown in Fig. [S] 
(straight bold lines on the upper and lower graphs; there 
also the dependencies for the case n — 10 are shown). Let 
us note that these dependencies (bold lines) also corre- 
spond to the case t = independently of the parameter k 
value (see also Eq. (jT5"]) ). 
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FIG. 6: (Color online) Dependencies of the group velocity and 
intensity of the transmitted light on the temperature for the 
two-level system with different values of the parameter k. 



Now we use the expressions that define the values of 
the group velocity and intensity of the transmitted light 



v s (t) = c{n'(t)+Lu[dn'(t)/du}Y 
I(t) = I exp[-n"(t)kL], 



(21) 



where L is the characteristic size of the atomic cloud (in 
calculations we set L = 0.004 cm). For the mentioned 
physical characteristics, one can obtain the dependen- 
cies corresponding to the signals tuned up exactly to the 
resonant frequency (A = 0). The corresponding curves 
are shown in Fig. [5] As one can conclude from Fig. [51 
the propagation velocity of the light pulse at the tem- 
peratures lower than the critical temperature in the case 
k«1 (as it is expected) practically does not depend on 
the temperature. 

Note that inequality k«1 usually takes place for di- 
lute gases of alkali-metal atoms that interact with a laser 
field (k « 0.025 for the parameters of the experiment [l[). 
Probably, this inequality takes place for the most BEC- 
related experiments at the present moment. Really, as 
it comes from the definition (fT2f of the parameter k, to 
increase its value in one order of magnitude, the density 
of atoms in a condensate must be increased by 3 orders 
of magnitude. This requirement results in the instabil- 
ity of the BEC phase due to an increase of the number 
of three-body collisions (see Ref. [12] )• Also it is easy to 
see that the use of the high-frequency radiation (with the 
large wave number k) cannot result to the considerable 
increase of the parameter k. It comes from the fact that 
the high-frequency transitions correspond to the high lev- 
els of the atomic spectrum. But the linewidth of these 
levels is much greater than the linewidth corresponding 
to the lower states. Hence, probably, the only way to 
a significant increase of the k value is to use the states 
with the smaller linewidth 7 (see Eq. (jT2"j) ). It means 
that to increase the mentioned value, one needs to select 



the levels whose probability of a spontaneous transition 
is much less than was used in the above calculations. 

Therefore, the cases K > 1 that are shown in Fig. [5] 
correspond to the "long-living" levels. One may relate 
these states to the levels with the forbidden dipole tran- 
sitions (e.g., hyperfine levels of the ground state, see in 
this case Ref. @) or the sublevels, whose relative transi- 
tion intensity is rather low. One can conclude that in this 
case, the temperature effects must have a strong impact 
on the ultraslow light phenomenon in a BEC. 



B. Temperature effects in three-level systems 

Note that two-level systems from the standpoint of 
the experiments dealing with the ultraslow-light phe- 
nomenon in a BEC can be inconvenient. As it is known, 
this phenomenon is realized mostly in three-level systems. 
These systems have some characteristic advantages: the 
absorption rate of the signal can be decreased by a special 
tuning of the laser frequency, one can get positive time 
delays of the pulses (positive sign of the group velocity) , 
and also one can use the magnetic field to control the 
group velocity of the ultraslow pulses [9( or use an addi- 
tional coupling laser to provide the electromagnetically 
induced transparency [1 31 ] - 

Therefore, in the framework of the developed approach 
let us study the temperature dependencies of the group 
velocity and intensity of the transmitted light in a three- 
level system. To this end, let us consider the system that 
is schematically illustrated in Fig. [7] (upper part). 

Note that in the framework of the developed approach, 
we do not account quantum interferences in the system. 
Therefore, in this case, the derived above expressions 
change insignificantly. In particular, due to the additive 
contribution of all quantum states to the Green functions 
(see Eq. ([6])), the functions G\ and G2 characterizing the 
response of all particles in the system can be written as 
follows (cf. Eq. (P]l): 



G 2 



F'(A-9,/2)] 
I'(A-Q/2,t;b)}, 



--(l-t 3 / 2 )[F'(A + n/2) 

+ bt[l'(A + Q/2,t;b) 
=(1 - t 3 / 2 )[F"{A + ft/2) + F"(A - Q/2)} 

+ bt[l"(A + Q/2, t; b) + I" (A ~ Q/2, t; b)], 



(22) 



where Q = Ae mag /j and Ae mag is the energy spacing 
between the excited states. Here and below, we con- 
sider that the splitted levels belong to the same multi- 
plet, i.e., |cr 12 | 2 /"7i2 = |cri3 1 2 /7i3- Hence, in accordance 
with Eq. (|20p . we can set 012 = 043 = a. Thus the for- 
mulas (Unj), (fTT)) . and (f2"Tj) that define the dependence 
of the dispersion characteristics on the temperature have 
the same form. 

To get the numerical estimates, let us consider the 
condensed gas of sodium atoms with the density v = 
5 x 10 12 cm -3 . For this gas, we get a = 0.052, n = 0.016, 
and b = 13.46. Also we consider that there are two 
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FIG. 7: (Color online) Scheme and dispersion dependencies 
of the three-level system in a BEC state. The used value for 
the parameter a — 0.052. 



nearby excited levels in the system (states |2) and |3) 
in Fig. [7J. We also set the spacing between them sev- 
eral times larger than the linewidth 7, Ae mag = 87. For 
this system, one can get the graphs characterizing the 
dispersion characteristics of the system that are shown 
in Fig. [7J (bold curves). Note that in this case (k -C 1), 
the functions n'(u>) and n"(u>) practically do not depend 
on the temperature. Let us emphasize that we regard 
the system as a superposition of two two-level systems 
(Autler-Townes treatment). This is the reason that the 
absorption rate does not go to zero at A = (Fig. [Jj bot- 
tom) as in the electromagnetically-induced transparency 
(EIT) regime P, [3. 

Within the developed approach, one can also study 
the dependencies in the case when the temperature ef- 
fects have a significant influence on the system response 
(as it is mentioned above, this corresponds to the case 
k > 1). In particular, taking n — 10 and the same val- 
ues of the other parameters, we get the dispersion curves 
at different temperatures that are shown in Fig. [7] It 
should be noted that in this case, the steepness of the 
slope of the refractive index decreases with the temper- 
ature. At the same time, the absorption in the central 
region increases with the temperature. This behavior can 
be explained by the fact that the density distribution of 
the non-condensate particles by the momentum p in the 
case ac > 1 is not so "sharp" as for the condensate par- 
ticles (see Eq. (JTJ). As a result, the absorption rate of 
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FIG. 8: (Color online) Dependencies of the group velocity and 
intensity of the transmitted light on the temperature for the 
three-level system with different values of the parameter k. 



the pulses tuned up to the resonance (A < 1) decreases 
and the absorption rate of the pulses detuned from the 
resonance (A > 1) increases (see Figs. 21 [5j and[7J with 
the temperature. 

The mentioned behavior of the main macroscopic pa- 
rameters with an account of Eq. ([2^]) results in the de- 
pendencies that are shown in Fig. [HI As one can see 
from the graphs, the velocity of a light pulse at the tem- 
perature lower than the critical temperature in the case 
k <C 1 (analogously to the two-level case) is practically 
constant. It corresponds to the fact that the refractive 
index profile and absorption rate in this case changes in- 
significantly. Let us emphasize that in this case for the 
calculations we used the same parameters of a gas with 
a BEC as were realized in Ref. jl| . 

At k > 1, the situation differs from the two-level case. 
Let us note that here one can see that both the group 
velocity and absorption rate of the pulse increase with the 
temperature. It can be explained by the fact that for the 
detuned pulses the steepness of the slope of the refractive 
index decreases while the absorption increases with the 
temperature (see Fig. [JJ. Thus the system becomes not 
so convenient for the realization of the ultraslow light 
phenomenon as in the zero-temperature limit. 



IV. CONCLUSION 

In this paper, we studied some optical properties of di- 
lute gases in the BEC state at finite temperatures. The 
analysis is based on the Green-function formalism. We 
found analytical expressions for the scalar Green func- 
tions that characterize a linear response of a gas to an 
external electromagnetic field (laser). We studied the 
characteristic dependencies of these functions on the tem- 
perature and frequency detuning. We made a comparison 
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of the contribution of the condensate and non-condensate 
particles into the effect. The ultraslow-light phenomenon 
in a BEC is studied both on the examples of two- level and 
three-level systems. 

In particular, it is shown that for a light pulse, which 
is tuned up close to the dipole-allowed transitions, the 
dispersion characteristics of a gas of alkali-metal atoms 
weakly depend on the temperature in the region from 
the absolute zero to the critical temperature. Therefore, 
the group velocity of the pulses in this system weakly 
depends on the temperature. This fact allowed us to 
conclude that the results of the previous works, where 
the authors studied the response of a gas in the limit of 
zero temperatures, T — > 0, are correct and can be used 
also in the region < T < T c . 

It is significant to note that in the experiment [l[ , data 
relating to the dependence of the group velocity of light 
pulses on the temperature of an ultracold gas of sodium 
atoms were obtained. From these results, one can con- 
clude that the group velocity weakly depends on the tem- 
perature in the region < T < T c . In the present paper, 
in the case of a three-level system, we use the same pa- 
rameters of a gas as in Ref. [lj and we come to the same 
conclusion in our research. But the statement that the 
results of our theory are verified by the mentioned ex- 
periment is not quite correct. It should be noted that 
in Ref. |l[ , the effect of the electromagnetically induced 
transparency was used. But in the framework of the in- 
troduced approach that is based on the Green-functions 
formalism, we cannot account for quantum interference 



effects in the system (see in this case Ref. [3(). Proba- 
bly, the similarity of the results may be a sign that the 
quantum interference effects in the mentioned experiment 
have a weak influence on the dependence of the group ve- 
locity on the temperature of a gas with a BEC. But, obvi- 
ously, this statement requires an additional experimental 
verification. 

We also studied the cases when the temperature effects 
can have a strong impact on the dispersion characteris- 
tics of gases with a BEC. In our opinion, this situation 
may be realized when the frequency of an external field 
is tuned close to the transitions between "long-living" 
states. These states correspond to the excited levels with 
the forbidden dipole transitions (upper hyperfine levels of 
the ground state of alkali-metal atoms) or the levels with 
a low relative intensity of the transition. In this case, 
as one can conclude from the results of the paper, the 
parameters of slowing impair with the temperature in- 
crease. Thus, the achievement of the lower temperatures 
of a gas is more necessary in this case. 
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